Lattice pinning of magnetic domains in the helimagnet Ba2CuGe2C>7 
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The layered magnetic compound Ba2CuGe2C>7 exhibits spiral antiferromagnetic order thanks to 
a Dzyaloshinskii-Moriya (DM) anisotropy that is allowed by crystal symmetry. Here we theoretically 
examine some finer issues such as the experimentally observed lattice pinning of the propagation 
f****) , vector of helical magnetic domains along the crystallographic (1,1,0) or (1,1,0) direction. We find 

that DM anisotropy alone would actually lead to pinning along the (1,0,0) or (0,1,0) direction, but 
£Nj ' agreement with experiment is restored upon including an additional exchange anisotropy that is also 

> | , consistent with symmetry. The present results also shed light on the so-called bisection rule which 

O j 1 has been abstracted from experiment in presence of an in-plane magnetic field. 
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PACS numbers: 75.30.Ds, 75.30.Gw 



Ba2CuGe2C>7 is an essentially two-dimensional spin system that exhibits spiral antiferromagnetic order due to a 
Dzyaloshinskii-Moriya (DM) anisotropy [1,2]. A series of experiments, including standard magnetometry, elastic, and 
inelastic neutron scattering [3-7], have established that the ground state is an incommensurate spiral with period 
^ (pitch) L w 37 lattice units (219 A) at zero external magnetic field. When a field of strength H is applied along the 
c-axis the period L = L(H) grows to infinity at H c w 2 T and the ground state degenerates into a commensurate spin- 
flop state for H > H c . Thus one obtained the first experimental realization of an incommensurate-to-commensurate 
(1C) phase transition whose gross features were theoretically predicted some time ago by Dzyaloshinskii [8]. 

Ba2CuGe207 is an insulator and its magnetic properties may be understood in terms of localized spins. Crystal 
symmetry provides the necessary input for the construction of a discrete Heisenberg Hamiltonian, but theoretical 
analysis is greatly facilitated by a continuum approximation which leads to a certain effective field theory that is a 
variant of the relativistic nonlinear a model. The continuum approximation is justified by the fact that anisotropy 
is sufficiently weak and thus leads to a ground state which is a spiral with a reasonably large period. Hence the 
nonlinear a model may safely be employed for the study of ground-state properties and the corresponding low-energy 
dynamics. The resulting theoretical picture was found to be in general agreement with experiment [3-7] and further 
' clarified the precise nature of the Dzyaloshinskii-type IC transition [9,10]. Yet, discreteness effects are important to 
, understanding some finer issues of experimental interest and are the main focus of the present paper. 
fSJ ' The unit cell of Ba2CuGe207 is partially illustrated in Fig. 1 where we display only the magnetic Cu 2+ ions with 
\ spin s = 1/2. The lattice constants are a = b = 8.466 A and c = 5.445 A. Since the Cu atoms form a perfect square 
lattice within each layer, it is also useful to consider the orthogonal axes x, y and z obtained from the conventional 
crystal axes a, b and c by a 45° azimuthal rotation. The major spin interaction between nearest in-plane neighbors 
along the x or y axis is antiferromagnetic, while the interaction between nearest out-of-plane neighbors along the z 
axis is ferromagnetic and rather weak [3]. Therefore, we are effectively dealing with a 2D antiferromagnet defined on 
. a square lattice with natural axes x and y and lattice constant d = a/\/2 = 5.986 A. 

Hence, to a first approximation, the spin dynamics is described in terms of a 2D isotropic Heisenberg model. 
However, because of the low tetragonal symmetry of this crystal (space group P42im), the effective Heisenberg 
Hamiltonian may contain an interesting combination of antisymmetric (DM) and symmetric exchange anisotropics. 
In the following, anisotropy is introduced in steps of increasing complexity, invoking the results of a complete symmetry 
analysis [9,10] as they become necessary. 

Experiments were initially [3-6] analyzed in terms of a relatively simple extension of the isotropic Heisenberg model 



defined by the Hamiltonian 

■ 

W=Y. [</(Sr-S r 0+D rr , -(SrXSrO], (1) 
£^ ' <rr'> 

c3 ■ 

where the sum extends over all bonds < rr' > connecting any two neighboring sites r and r' on a square lattice. 
Symmetry requires that the exchange constant J is the same for all bonds, while the constant vectors D rr < which 
account for pure DM anisotropy are of the form 

D rr ' = (0, D±,±D Z ) for bonds along x, 

D rr , = (D ± ,0,±D z ) for bonds along y, (2) 

where Cartesian components are taken along the axes x, y and z of Fig. 1, and a sign alternation at consecutive 
bonds is present in the third component ±D Z . In fact, all available experimental evidence suggests that -D Z <C-D± in 
Ba 2 CuGe 2 07. 
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FIG. 1. Partial illustration of the unit cell of Ba2CuGe207 displaying only the magnetic Cu 2+ ions, denoted by solid circles, 
which form a perfect square lattice within each layer. 

Therefore, while possible implications of a finite D z are studied in Refs. [9,10], we set D z = throughout the 
present paper. Hamiltonian (1) is then written in a completely explicit form: 

W = ^2 [ J S m . n ■ ( S TO+ i,„ + S mi „+i) + Dj_(S^ n n S^ n+l n — S^ i n S^ n+l n ) 

mn 

+ D±_{ S^ n n S m n+1 — S m n S^ n n+1 ) , (3) 

where m and n are integers that advance along the x and y axes of the square lattice and S m-n = (S^ n , n , n ) 
is the spin operator at site (m, n). At this point we invoke the classical approximation whereby spin operators are 
treated as classical vectors with magnitude „ = s 2 . The ground state is then described by the classical spin 
configuration that minimizes Hamiltonian (3). To search for the minimum we employ the spherical parametrization 



Q X 


= *(- 


-|\m+7i 


sin6 mi „ 


cos $ m ,„ 


cy 

^ m , n 


= *(- 


-|\m+7i 


sin6 mi „ 


sin $ m ,„ 


QZ 


= *(- 




COS ®rn,n 


1 



(4) 

where the staggering factor (— l) m +™ reflects the antiferromagnetic nature of the basic exchange interaction. Inserting 
Eq. (4) in Eq. (3) yields an energy functional whose variation leads to a coupled system of nonlinear equations for 
the angular variables Q m , n and <fr m ,n- We shall not write out these equations explicitly, but merely state a class of 
analytical spiral-like solutions: 

6 m ,™ = m 7i +7172, $ m , n = -V>, (5) 

where the constants 71, 72, and ip are constrained only by the curious relation 

tan^^L. (6) 
sin 71 



w = s J 



The corresponding energy per unit cell is given by 

2 — cos 71 — cos 72 — e (sin 71 cos ip + sin 72 sin ip) (7) 

where the trivial constant 2s 2 J amounts to removing the energy of the pure Neel state, and 

e = DJJ (8) 
is a dimensionless anisotropy constant whose magnitude is left arbitrary for the moment. 
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The helical spin configurations constructed above are stationary points of the energy functional for any choice of 
the constants 71, 72, and ip that satisfy the constraint (6). However, the absolute minimum of the energy is achieved 
for the specific choice 71 = 72 = arctan (e/y/2) and ip = 7r /4, or 



(m + n) arctan (e/\/2), 



1>, 



-tt/4, 



w 



2s 2 J 



1- v/l + e 2 /2 



2 2 7 
= e s J 



e 

16 



(9) 



which is a screw-type spiral whose propagation vector points along the crystalographic b = (0,1,0) axis, while the 
spin rotates in a plane perpendicular to the same axis. An equivalent solution that propagates along the a = (1, 0, 0) 
axis is also possible but will not be given further consideration. 

The preceding result is already contradicted by the experimental fact that the spiral propagates along the x = 
(1,1,0) or the y = (1,1,0) axis [3-7]. To make this fact completely explicit we consider a special member of the 
family of solutions (5) obtained by setting ip = 0, and thus 72 = by virtue of the constraint (6), while the remaining 
constant 71 is chosen to minimize the reduced energy w = s 2 J(l — cos 71 — esin7i) to yield 71 = arctan (e) and 



m arctan (e), 



0. 



s 2 J 



e 2 s 2 J 



1 e z 

-2 + ¥ + - 



(10) 



which is a spiral that propagates along the x axis while the spin rotates within the xz plane. Both of these facts are 
consistent with experiment, but the energy of the spiral of Eq. (10) is larger than the energy of the screw-type spiral 
of of Eq. (9) to within terms of order e 4 . In other words, the pure DM Hamiltonian (3) cannot explain the observed 
pinning of the propagation vector of the spiral along the x or y axis. 

Before discussing a resolution of this apparent discrepancy, it is instructive to view the analytical solutions given by 
Eqs. (5)-(7) from a slightly different perspective. Thus we attempt to first minimize the energy (7) at fixed angle ip, 
with 71 and 72 related by the constraint (6). The result of this minimization is written here in a perturbative form: 



sin 71 = e cos ip 
sin 72 = e sin ip 



1 — - (cos 4 V> + sin 4 ?/;) e 2 + ... 
1 — i (cos 4 ?/; + sin 4 ?/;) e 2 + ... 



and the corresponding energy per unit cell is given by 

w = £ 2 s 2 J [-\ + ^ [3 + cos(4V0]e 2 + ...}. 



(11) 



(12) 



This result is consistent with Eqs. (9) and (10), for ip = 7r/4 and ip = 0, respectively, and reinforces our earlier 
conclusion that the absolute minimum of the energy is achieved for ip = tt/4 which leads to the screw- type spiral of 
Eq. (9). 

For sufficiently weak anisotropy (e<§Cl) Eq. (11) may be further approximated by 71 w £cosip and 72 ~ £simp 
which are inserted in Eq.(5) to yield 



£ (m cos ip + n sin ip) , 



= -ip. 



(13) 



This is a class of degenerate helical configurations with common energy w w — £ s J/2 which is independent of angle 
ip. Eq. (13) describes a helix that propagates along the x' axis of Fig. 2, which forms an angle ip with respect to 
the x axis, while the spin rotates in the x" z plane whose intersection with the basal plane forms an angle — ip with 
respect to the x axis. The angle formed by the normal to the spin plane (axis y") and the spiral propagation vector 
(axis x') is bisected by the conventional crystal axis b = (0, 1,0). Therefore, to leading e -C 1 approximation, theory 
predicts that spirals may propagate in an arbitrary direction (angle ip) as long as they conform with the "bisection 
rule" just described. 

However, when the terms of order e 4 are included in the energy, as in Eq. (12), a screw-type spiral propagating 
along the a or b axis (ip — ±7r/4) is predicted to be energetically favorable. As mentioned already, this prediction 
is contradicted by experiment where spirals always propagate along the x or y axis (ip — or n/2) in the absence 
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of an in-plane external magnetic field. Therefore, the pure DM anisotropy included in Hamiltonian (1) needs to be 
ammendcd. 

A more general type of anisotropy is suggested by the early work of Kaplan [11], and of Shckhtman, Entin-Wohlman, 
and Aharony [12-14], which is usually referred to as the KSEA anisotropy. The latter was already invoked in Refs. 
[7,9,10] to explain the observed distortion of the spiral from its ideal sinusoidal shape, as well as the structure of the 
magnon spectrum determined via inelastic neutron scattering. In the continuation of this paper we show that the 
KSEA anisotropy also explains: (a) the observed lattice pinning of the spiral propagation vector along the x = (1, 1, 0) 
or y — (1,1,0) direction in the absence of an in-plane field, and (b) a remnant of the bisection rule observed via neutron 
diffraction in the presence of an in-plane field [5] . The KSEA extension of Hamiltonian (3) reads 



W 



where 



and 



E[ j qx OI l7Cy CV I 7 q z oz 

mn 

I J qx qx I T qy qy _i_ J qz qz 

' J 2 '-'iTLti^m.ji+l "r" J l D m,n D m,n+l J 3 J m ,n°m,n+1 

_\ J~) . ( q z q x _ q x q z i qv q z _ q z qV ) (ia) 

Ji = J - A, J 2 = J + A, J 3 = J - A (15) 



A = ff = \*J. (!6) 

One should keep in mind that Hamiltonian (14) is still consistent with the underlying space group P\2\m but is not 
the most general Hamiltonian allowed by symmetry [9,10]. Nevertheless, Eq. (14) has been the starting point for a 
reasonably successful analysis of a vast set of experimental data and will be adopted here without further questioning. 

Our aim is then to find the classical minima of Hamiltonian (14). One may again employ the spherical parametriza- 
tion (4) but the resulting coupled equations for O mi „ and $ mi „ do not seem to admit an analytical solution (when 
A 7^ 0) of the type obtained earlier in Eqs. (5)-(7) within the pure DM model (A — 0). Instead, our strategy is to 
derive a perturbative expansion for small values of e = D±/J of practical interest. 

To leading order in e, the ground state may be obtained by minimizing a suitable continuum energy functional 
[7,9,10]. The resulting spiral generalizes Eq. (13) to the extent that anharmonic distortions occur in the profile of the 
angular variable Q m ,n which is no longer a linear function of e(mcos?/' + nsinip) s=y (a; cost/; + ys'mtp). Instead, the 
complete solution may be written as 

©m,« ~ (x cos V> + y sin tp) , $ m ,„ = -tp, (17) 

where tp is an arbitrary constant angle and 9(u) is a function of u = x cos ip + ysintp that may be expressed in terms 
of an elliptic integral. This family of solutions is degenerate in the sense that the energy is independent of the angle 
ip. As a consequence, the bisection rule summarized in Fig. 2 remains valid in the presence of KSEA anisotropy to 
leading order in e. 

However, the above degeneracy is again broken beyond the leading approximation. To demonstrate this fact we 
have systematically carried out the continuum expansion beyond the leading approximation. Technical details are too 
cumbersome to be described here in any detail, but the main conclusion may be stated by quoting the final result for 
the ground-state energy per unit cell calculated to order s 4 : 

w = e 2 s 2 j{ - V + C* + C\ cos (4V0 s 2 + ...}. (18) 

When this procedure is applied to the pure DM model (A = 0) we find that <5 2 = 1, Co — 3/32, and C\ = 1/32, 
which agree with the results of Eq. (12) and thus provide an important check of consistency. In the presence of KSEA 
anisotropy (A = e 2 J/4) we find the values 

5 2 = 0.53189772, C = 0.00098724, C\ = -0.00474723, (19) 

obtained by numerical calculation of a large number of elliptic-type integrals. The most important feature of Eq. (19) 
is that the constant Cl is now negative and thus leads to minimum energy in Eq. (18) for ip = or 7r/2; i.e., for a 
spiral that propagates along the x or y axis, in agreement with experiment. 



4 




FIG. 2. Illustration of the ideal bisection rule. The spiral propagation vector points along the x' axis, while the spin rotates 
within the x"z plane. This ideal rule is experimentally realized when a sufficiently strong in-plane magnetic field H is applied 
along the y" axis. 



The strength of anisotropy appropriate for the description of Ba2CuGe20y is estimated from the value of the spiral 
pitch measured at zero field [7] to be e = D±/J s=s 0.18. Therefore, anisotropy appears to be sufficiently weak for the 
validity of the leading-order continuum approximation. For instance, the relative correction induced by the second- 
order term in Eq. (18) is of the order 5 x 10~ 4 . Nevertheless, such a tiny correction is important for explaining the 
correct lattice pinning of the spiral propagation vector along the x or y axis, as well as some measurable deviations 
from the ideal bisection rule summarized in Fig. 2. 

The bisection rule was initially discovered through neutron diffraction [5] in the presence of an in-plane magnetic 
field 

H = (cos X ,sinx,0)ff, (20) 

which forms an angle % with respect to the x axis. The spiral was then observed to adjust its propagation vector 
according to the rule 

X + V> = |, (21) 

in order to minimize the Zeeman energy induced by the magnetic field. Returning to Fig. 2, a spiral that initially 
propagates along the x axis (ip = 0) is reoriented to propagate along the x 1 axis (ip ^ 0) so that the normal to the 
spin plane (axis y") points along the applied magnetic field. 

Actually, the above rule was observed only for sufficiently strong fields H > 0.5 T. For weak fields, reorientation of 
the spin structure to conform with the ideal bisection rule is inhibited by the small tetragonal anisotropy that is present 
in the second order term of Eq. (18). The appropriate modification of Eq. (21) is then derived by an argument similar 
to that given in Ref. [5], now taking into account the theoretical prediction for the tetragonal anisotropy contained 
in Eq. (18): 

g2= 8Asin W ) 

sin2(^ + X )' V ' 



where 



2sJ V 2 Idle 4 
9ab^Bj (sin 2 6>) 



(23) 



is a numerical factor which is here calculated in terms of the constant C\ of Eq. (19), the average (sin 2 6>) = 0.56517 
taken over a period of the spiral profile of Eq. (17), and some basic parameters determined by independent experiments 



5 



[3-7]; namely, the spin value s = 1/2, the exchange constant J = 0.96 meV, the DM anisotropy e = D±/J = 0.1774, 
the gyromagnetic ratio g a b = 2.044 for a field applied in some direction within the basal plane, and the Bohr magneton 
Ms- 




FIG. 3. Deviations from the ideal bisection rule predicted by Eq.(22) applied for A = 1.1 x 10~ 3 T 2 (dotted lines) and 
A — 1.95 x 1CP 3 T 2 (solid lines). Note that the ideal bisection rule (x + ip = 7r/2) is indeed realized for a sufficiently strong 
in-plane magnetic field H > 0.5 T. Experimental data denoted by symbols were extracted form Fig. 3 of Ref. [5]. 

Putting everything together, we find that A = 1.1 x 10~ 3 T 2 . The corresponding predictions of Eq. (22) are 
depicted by dotted lines in Fig. 3 where we display the direction of the spiral propagation (angle ip) as a function of 
applied field (H) for various field orientations (angle x). The ideal bisection rule of Eq. (21) is clearly reproduced 
for sufficiently strong fields (H > 0.5 T) but measurable deviations occur for weak fields. Altogether, our theoretical 
predictions are in reasonable agreement with experiment. A better fit is obtained by treating the constant A in 
Eq. (22) as an adjustable phenomenological parameter determined by a global fit of the experimental data, as was 
originally done in Ref. [5] to obtain A = 1.95 x 10~ 3 T 2 . The corresponding predictions of Eq. (22) are also depicted 
in Fig. 3 (solid lines) and provide an improved overall fit of the data. 

To summarize, the KSEA extension of the pure DM anisotropy is capable of explaining a large set of data pertaining 
to the IC phase transition observed in Ba2CuGe207, as well as some finer issues such as the lattice pinning of helical 
magnetic domains. One should keep in mind that DM+KSEA is not the most general type of anisotropy allowed by 
symmetry [7,9,10]. In this respect, we mention here that a significant departure from KSEA seems to occur in the 
layered antifcrromagnet K 2 V30s. The space group of this compound is different (PAbm) but leads to a Heisenberg 
Hamiltonian that is very similar to the one encountered in Ba 2 CuGe207. However, an experimental investigation [15] 
of K2V3O8 revealed the occurrence of interesting spin- flop and spin- reorientation transitions but provided no evidence 
for incommensurate magnetism. A preliminary theoretical analysis [15,16] concluded that K2V3O8 is characterized 
by an easy-axis anisotropy, which is impossible to occur in the KSEA limit [10]. 

This paper is dedicated to our friend Professor Victor G. Bar'yakhtar on the occasion of his 75th birthday - 
helimagnetism has been one of the many subjects of interest during his long and fruitful career [17]. 
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